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ABSTRACT 

This  report  describes  CHAP,  a computer  language  for  processing  line 
drawings  represented  in  terms  of  the  8-direction  chain  code.  The  lan- 
guage consists  of  a set  of  FORTRAN  subroutines,  most  of  which  are  machine 
independent.  The  subroutines  may  be  called  by  a user  program  to  accom- 
lish  specific  line-drawing  processing  tasks.  The  chain  data  is  stored 
in  packed  (machine-dependent)  form  and  is  unpacked  for  processing.  The 
report  lists  all  currently  available  CHAP  subroutines  and  provides  all 
necessary  information  for  a CHAP  programmer.  Specific  information  for 
different  machine  implementations  are  given  in  the  appendices.  Thus  far 
an  implementation  exists  only  for  an  IBM  360  computer.  However, 
implementations  for  UNIVAC  1108/1110,  CDC  6600,  and  ADAGE  AGT  30/130  are 
presently  being  planned. 
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1.  INTRODUCTION 


Since  the  introduction  of  the  chain  code  representation  for  free- 
form line  drawings  in  1961  [l],  many  algorithms  for  processing  chain- 
encoded  line  drawings  have  been  developed  [2-4],  CHAP  is  a result  of  an 
effort  to  consolidate  these  algorithms  into  a single  line-drawing  pro- 
cessing language.  The  language  consists  of  a set  of  FORTRAN  subroutines 
which  can  be  invoked  by  users'  programs.  The  design  of  CHAP  is  analogous 
to  that  of  the  PAX  image  processing  language  [5].  The  CHAP  subroutines 
are  mainly  written  in  standard  ANSI  FORTRAN  IV  and  implemented  on  an  I3M 
360/67  computer.  A study  was  made  to  determine  how  CHAP  could  be  made 
readily  portable  to  other  computers.  At  issue  here  were  questions  con- 
cerning differences  in  word  length  (ranging  from  l6  to  60  bits),  subtle 
difference  in  the  particular  FORTRAN  versions,  variations  in  I/O  handling, 
and  problems  of  linking  to  other  image  processing  software. 
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2.  THE  CHAP  LANGUAGE 


A chain  is  a sequence  of  octal  digits  a^ , i = 1 n,  representing 
a curve  in  terms  of  a string  of  short,  connected,  straight-line  segments. 

Each  octal  digit  corresponds  to  a particular  directed  line  segment 
(chain  link)  of  length  T(/2)S  oriented  at  an  angle  of  m(^5°)  relative  to 
the  positive  x axis;  n is  the  octal  value  of  a.,  s is  the  modulo-two 
value  of  m,  and  T is  a scale  factor.  A chain  may  be  regarded  as  a sequence 
of  line  segments  connecting  nodes  that  are  8-adjacent  to  each  other  on  a 
square  lattice  of  spacing  T.  Even-valued  chain  links  connect  nodes  lying 
vertically  or  horizontally  adjacent  and  are  of  length  T;  odd-valued  links 
connect  nodes  lying  diagonally  adjacent  and  are  of  length  T/2.  The  scale 
factor  T is  normally  set  to  unity.  A chain  usually  contains  special 
signal  codes  (i.e.,  control  codes)  distinguished  by  the  link  sequence 
OULjLj.  Page  L.l  lists  the  presently  assigned  chain  signal  codes.1* 

In  CHAP,  chain  links  are  stored  in  two  formats — packed  and  unpacked. 

In  the  packed  format,  the  chain  links  are  stored  in  an  array  with  | n./3  ■ 
links  in  one  word,  where  n is  the  number  of  bits  per  word.  Thus  for  the 
CDC-6600  implementation,  twenty  links  would  be  stored  in  one  60-bit  word. 

A chain  begins  with  a chain  identification  number  (which  may  be  up  to  63 
octal  digits  long)  followed  by  the  x and  y coordinates  of  its  initium. 

For  example,  the  sequence  ' 0L07  0 5 100000L26000350U2700025 ' indicates 

that  chain  number  1000a  follows,  and  that  the  coordinates  of  its  initium 

are  x“358and  y= 250  A chain  is  terminated  by  the  end-of-chain  indicator  '0^00'. 

* A reader  unfamiliar  with  the  chain  coding  scheme  is  advised  to  read 
reference  U,  which  provides  a good  introductory  survey  of  the  subject. 
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CHAP  maintains  a catalog  of  all  chains  currently  in  memory.  The 
catalog  entry  for  a chain  contains  the  chain  identification  number  and  the 
pointers  to  the  first  and  last  memory  locations  of  the  chain.  This 
arrangement  allows  a user  to  specify  which  chains  are  to  be  processed  by 
the  various  routines  by  simply  referring  to  the  chain  identification 
numbers.  When  a chain  is  to  be  processed  by  a routine,  it  is  unpacked 
from  the  packed  format  into  a working  buffer.  Usually  only  a part  of  the 
chain  in  unpacked  at  a time,  processed,  and  then  the  next  part  of  the  chain 
is  unpacked.  There  are  three  working  buffers  in  CHAP.  Their  contents  are 
also  catalogued  by  CHAP.  This  avoids  unpacking  a chain  which  has  already 
been  unpacked  into  a working  buffer  by  a previous  call  to  a processing 
routine.  Since  some  routines,  such  as  chain  crosscorrelation  and  chain 
intersection,  operate  on  two  chains  simultaneously,  they  use  two  of  the 
working  buffers  at  the  same  time.  The  unpacking  routine  of  CHAP  stores  one 
link  per  word,  except  that  the  signal  codes  are  stored  with  the  two-digit 
code  LjL2  biased  V lOjowhen  the  signal  code  contains  no  data,  and  biased 
tylOOiowhen  it  does,  in  which  case  the  data  is  stored  in  the  next  memory 
locaticn(s).  Thus  the  chain  sequence  given  in  the  previous  paragraph  is 
unpacked  into  individual  words  as  follows  (Now  shown  as  FORTRAN  decimal  integers) 

107  chain-identification-number  code  (lOOio  +07) 

1*096  chain  identification  number  (I0000e=l*96io ) 

122  x-coordinate  code  (lOOi  o+26s=122i 0 ) 

29  x-coordinate  (35s=29io) 

123  y-coordinate  code  (lOOi  o+27s=123i 0 ) 

21  y-coordinate  (25»=21io) 

This  format  of  ’unpacked  working  buffers  allows  efficient  processing 
because  values  from  0 to  7 indicate  chain  directional  links  and  all  other 
values  indicate  signal  codes;  those  greater  than  lOOi  0 also  nave  one  or  more 
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following  data  words.  Note  that  links  preceded  by  a repeat  signal  code 
(such  as  0Ul7 ) are  completely  unpacked  into  the  number  of  links  indicated 
by  the  repeat  code  and,  therefore,  the  signal  code  does  not  appear  in  the 
unpacked  form. 

CHAP  can  also  be  used  to  create  a chain.  First  a working  buffer 
must  be  initialized  and  a chain  identification  number  must  be  assigned 
and  entered  into  the  buffer.  This  is  accomplished  by  the  CHAP  subroutine 
BUILD.  Then  links  or  sequences  of  links  are  entered  and  appended  to  the 
identification  number  and  any  signal  codes  in  the  buffer.  The  chain  may 
then  be  put  into  packed  form  and  catalogued  as  a CHAP  chain,  or  it  may  be 
processed  by  any  of  the  CHAP  chain  processing  routines. 
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3.  DEFINITIONS  OF  CHAP  TERMS 


Binary  Plane  - A two-dimensional  array  of  bits  which  are  set  to  1 at 
locations  of  chain  vertex  points  and  0 everywhere  else.  The  CHAP  binary 
plane  is  stored  in  common  block/PLANES/NX,NY,PLANE(NX*NY/NBW)  where  NX  is 
the  number  of  bits  (points)  in  the  x direction  of  the  plane  and  NY  is  the 

number  of  bits  (points)  in  the  y direction  of  the  plane.  PLANE  is  a 

linear  array  in  which  the  bits  of  the  plane  are  stored.  The  dimensic-  of 

PLANE  is  [NX*HY/NBW]  where  NBW  is  the  number  of  bits  per  word. 

Chain  - A sequence  of  octal  digits  called  links.  A chain  consists  of 

signal  codes  and  directional  links.  A chain  begins  with  a signal  code  giving  its 

identification  number,  usually  followed  by  signal  codes  giving  the  initial 

x and  y coordinates  of  the  chain.  The  initial  coordinates  are  set  zo  zero  if  they 

are  not  specified.  A chain  must  be  terminated  by  the  end-of-chain  signal  code  'OUOO'. 

Chain  data  file  - A data  file  containing  chains  in  the  following  format: 


olumn 

1-5 

Chain  identification  number 

6 

Ignored 

7-72 

Chain  data 

73-end 

Sequence  numbers 

A data  file  is  terminated  by  a zero  chain  identification  number. 

The  chain  data  file  may  contain  any  characters  but  only  octal 
digits  (0  through  7)  are  read  as  chain  links.  This  allows  users  to  freely 
insert  blanks  or  other  (non-octal-digit)  characters  to  improve  the  read- 
ability of  chain  data  files. 


Chain  identification  number  - A positive  integer  by  which  a chain 
is  uniquely  identified  to  CHAP.  In  chain  data  files,  it  is  a five-digit 
decimal  number  entered  in  columns  1-5.  It  may  be  assigned  by  the  user 
when  building  a chain  in  CHAP. 

Signal  codes  - In  a chain,  a sequence  of  links  beginning  with  OU  that 
indicates  the  start  of  control  information.  A list  of  all  currently 
assigned  signal  codes  is  given  in  Section  U. 

Directional  link  - A directed  straight-line  segment  of  length  1 or  72 
connecting  the  current  node  in  a square  lattice  of  unit  spacing  to  one  of 
its  eight  neighboring  nodes.  A link  is  identified  according  to  its  direc- 
tion by  a number  from  0 to  7. 


Link  - An  octal  digit  within  a chain. 

Packed  chain  - A chain  whose  links  are  packed  into  memory  words  with  l_n /3 J 
links  per  word,  where  n is  the  number  of  bits  per  word  of  the  particular 
computer  being  used. 

Unpacked  chain  - A chain  whose  directional  links  are  stored  one  per  word. 
Signal  codes  that  do  not  contain  any  data  are  biased  by  10jo  to  differentiate 
them  from  directional  links.  Signal  codes  that  contain  data  are  biased  by 
lOOio  to  differentiate  then  from  directional  links  and  signal  codes  without 
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data.  The  data  is  stored  in  the  locations  following  the  signal  code.  For 

example,  the  chain  ' 0^070501000  0h2600010  0^2700021  1251*  0^00'  is  unpacked  into 

individual  words  as  follows  (Now  shown  as  FORTRAN  decimal  integers): 

107  chain-identification-number  signal  code  (07+100io) 

512  chain  identification  number  (I000e=512i o ) 

122  x-coordinate  signal  code  (26  +100! 0=122i 0 ) 

8 x-coordinate  (108=8io) 

123  y-coordinate  signal  code  (278+10010=123i 0) 

17  y-coordinate  (2le=17io) 

1 directional  link 

2 directional  link 

5 directional  link 

« directional  link 

10  end-of -chain  signal  code  (00+1010) 


Chain  Vertex  - A chain  vertex  is  a gird-node  point  through  which  the  chain 
passes.  It  is  normally  denoted  by  the  number  of  the  link  which  departs 
from  it.  Thus  vertex  1 is  the  initium  of  the  chain,  and  the  last  link 
joins  vertex  n to  vertex  n+1.  Vertex  n+1  is  the  terminus  of  the  chain. 
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U.  LIST  OF  CHAP  CHAIN  SIGNAL  CODES 


Code 

ouoo 

0U01 

0h02 

0U03 

0^04 
01*0  5XYZ 
OUOTUV 
0U10WXYZ 
01*11 

01*127777 

JU13UWXYZ 

OUlUUVWXYZ 

0U15UVWXYZS 

Ol*17UXYZ 

0U20UN 

0l*21TUVWXYZ 

0l*22U 

0U23WXYZ 

0l*2l*XYZ 

0U25U 

0l*26VWXYZ 

Ol*27’AOCYZ 


Description 

End-o f-chain  indicator. 

Invisible  chain  follows. 

Visible  chain  follows.  Every  chain  is  assumed 
to  be  visible  initially. 

Remainder  of  word  is  to  be  skipped.  This  code 
is  ignored  by  CHAP. 

Valid  0U  sequence  of  directional  links 
Point  marker  number  (XYZ)a 

Chain  identification  number  (UV)8  digits  long  follows. 

Non-chain  data(WXYZ)8  digits  long  follows. 

Non-chain  data  follows  until  terminated  by 
end-of-ncn-chain-data  indicator . 

End-of-non-chain-data  indicator. 

Chain  node  number  (WXYZ)e  connects  with  (U)8  other  chains. 

Rotation  indicator. 

Scale  change  indicator. 

Link  U is  to  be  repeated  (XYZ)s  times. 

Link  U is  to  be  repeated  as  many  times  as  specified 
by  the  next  (N)g  ♦ 1*  digits. 

The  sequence  consisting  of  the  immediately  following 
(TUV)g  links  is  tc  be  repeat'd  (WXYZ)g  times. 

Color  indicator. 

Elevation  value  indicator  (WXYZ)8 
Grey-level  value  indicator  (XYZ)8 
Check  code.  U8  yields  modulo-8  sum  of  zero. 

Initial  x coordinate  is  (VWXYZ)g 
Initial  y coordinate  is  (VWXYZ)8 


l*.l 


5.  LIST  OF  SUBROUTINES  GROUPED  BY  FUNCTION 


CHAP  Initialization: 

CHAP3D  Block  data  program.  Allocate  and  initialize  the  CHAP 

labelled  common  blocks. 

Subroutines  used  for  chain  synthesis: 


BUILD 

XCOORD 

YCOORD 

LINK 

LINKSQ 

INVI3 

VISIBL 

CHLINE 

EOCH 


Start  to  build  a chain  by  initializing  buffer  and  inserting 
chain's  identification  number. 

Append  an  x coordinate  to  a chain. 

Append  a y coordinate  to  a chain. 

Append  a link  to  a chain. 

Append  a sequence  of  links  to  a chain. 

Append  invisible-chain- follows  signal  code. 

Append  visible-chain-follows  signal  code. 

Generate  a sequence  of  links  approximating  a straight  line. 
Append  end-of-chain  signal  code. 


Chain  management  subroutines: 


INPUT  Input  a chain  data  file. 

CHLIST  List  all  chains  in  CHAP  catalog. 

DELETE  Delete  a chain  from  CHAP. 

OUTPUT  Outputs  chains  to  I/O  unit  7 in  format  read  by  INPUT. 

ARRAY  Convert  a chain  into  an  array  containing  the  links  one  per  word. 

VERTEX  Convert  a chain  into  lists  of  vertex  coordinates. 

NPACK  Pack  a chain. 

UNPACK  Unpack  a chain. 

3PLANE  Clear  binary  plane  and  convert  a chain  to  binary  plane 

representation . 

APLANE  Convert  and  append  a chain  into  a binary  plane. 

BPRINT  Print  the  binary  plane. 

PRINT  Convert  chain  to  binary  plane  representation  and  print  the 

binary  plane. 

CHPAX  Convert  a CHAP  chain  to  a PAX  chain. 

CHPLOT  Plot  a chain  on  a digital  plotter. 


Link  management  routines: 

GTLINK  Unpack  a link  from  a packed  format  chain. 

OCTNUM  Form  an  octal  number  from  a sequence  of  links. 

STLINK  Store  a link  in  a packed  format  buffer. 
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Chain  analysis  subroutines: 


ANGLE 

AUTO 

CENTRD 

COARSE 

CROSS 

ECAREA 

INITXY 

INTERS 

LENGTH 

LMOM1 

LM0M2 

LNAREA 

LNCED 

MATCH 

MAXMIN 


MOM1 
M0M2 
PD  I ST 
POLYGN 

PNTCHD 

RESID 

ROSCAL 

WHEX 

XPROFL 

YPROFL 


Compute  angle  between  the  x axis  and  line  defined  by  two 
chain 

Compute  the  autocorrelation  of  a chain. 

Compute  the  x and  y coordinates  of  the  centroid  of  a chain. 
Requantize  chain  on  a coarser  grid. 

Compute  the  crosscorrelation  of  two  chains. 

Compute  the  enclosed  are  of  a closed  chain. 

Find  the  x and  y coordinates  of  the  chain  initium 
Find  the  intersection  of  two  chains. 

Compute  the  length  of  a chain. 

Compute  the  first  moment  of  a chain  about  a line. 

Compute  the  second  moment  of  a chain  about  a line. 

Compute  the  area  between  a chain  and  a line. 

Compute  the  max  and  min  distances  from  a chain  to  a line. 
Match  two  chain  profiles. 

Compute  the  max  and  min  extents  of  a chain  in  x and  y 
directions . 

Compute  the  first  moment  of  a chain  about  0,1*5,90,135  axes. 
Compute  the  second  moment  of  a chain  about  0,1*5,90,135  axes. 
Compute  the  distance  between  two  points  on  a chain. 

Determine  a polygonal  approximation  of  a chain  to  a specified 
precision. 

Compute  max  and  min  distances  from  a chain  to  a point . 
Compute  chain  residue. 

Rotate  and  scale  a chain. 

Compute  chain  width,  height,  1*5  and  135  degrees  extents. 
Compute  chain  profile  along  x axis. 

Compute  chain  profile  along  y axis. 
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ALPHABETICAL  LIST  OF  GLOBAL  NAMES  AND  COMMON  BLOCKS 


The  following  is  a list  of  all  global  nam.es  defined  in  CHAP  including 
Subroutines  (S),  Entry  to  a subroutine  (E),  Integer  functions  (IF),  and 
Block  data  C BD ) . 


Name  Tyre 


ANGLE  S 

APLANE  E 

ARRAY  S 

AUTO  S 

BPLANE  S 

BPRINT  E 

BUILD  S 

CENTRD  S 

CHAPBD  BD 

CHLINE  3 

CHLIST  S 

CHPAX  S 

CHPLOT  3 

COARSE  3 

CROSS  S 

DELETE  S 

ECAREA  S 

EOCH  S 

GTLINK  IF 

INITXY  3 

INPUT  S 

INTERS  S 

INVIS  S 

LENGTH  S 

LINK  , S 
LINKSQ  3 

LMOM1  S 

LM0M2  S 

LNAREA  S 

LNCHD  S 

MATCH  S 

MAXMIN  3 

MOM1  3 

M0M1A  E 

M0M2  S 

M0M2A  E 

NPACK  S 

OCTNUM  IF 

OUTPUT  S 

PDIST  S 


Description 


Compute  angle  between  the  x axis  and  the  line  defined  by 
two  chain  vertices. 

Convert  and  append  a chain  into  a binary  plane. 

Convert  chain  into  array  containing  the  links  one  per  word. 
Compute  the  autocorrelation  of  a chain. 

Clear  binary  plane  and  convert  a chain  into  it . 

Print  binary  plane. 

Start  to  build  a chain  by  initializing  buffer  and  inserting 
chain's  identification  number. 

Compute  the  x and  y coordinates  of  the  centroid  of  a chain. 
Allocate  and  initialize  the  CHAP  labelled  common  blocks. 
Generate  a sequence  of  links  approximating  a straight  line. 
List  all  chains  in  CHAP  catalog. 

Convert  a CHAP  chain  to  a PAX  chain. 

Plot  a chain  on  a digital  plotter. 

Requantize  chain  on  a coarser  grid. 

Compute  the  crosscorrelation  of  two  chains. 

Delete  a chain  from  CHAP. 

Compute  the  enclosed  area  of  a closed  chain. 

Append  end-of-chain  signal  code  to  a chain. 

Unpack  a link  from  a packed  format  chain. 

Find  the  x and  y coordinates  of  the  chain  initium 
Input  a chain  data  file. 

Find  the  intersections  of  two  chains. 

Append  invisible-chain-follows  signal  code. 

Compute  the  length  of  a chain. 

Append  a link  to  a chain. 

Append  a sequence  of  links  to  a chain. 

Compute  the  first  moment  of  a chain  about  a line. 

Compute  the  second  moment  of  a chain  about  a line. 

Compute  the  area  between  a chain  and  a line. 

Compute  the  max  and  min  distance  from  a chain  to  a line. 
Match  two  chain  profiles. 

Compute  the  max  and  min  extents  of  a chain  in  x and  y 
directions . 

Compute  the  first  of  a chain  about  0,^5,90,135  axes. 

Entry  to  routine  M0M2  (CHAP  internal  use  only). 

Compute  the  second  moment  of  a chain  about  0,^5,90,135  axes. 
Entry  to  routine  M0M1  (CHAP  internal  use  only). 

Pack  a chain. 

Form  an  octal  number  from  a sequence  of  links. 

Output  chains  to  I/O  unit  7 in  format  read  by  INPUT. 

Compute  the  distance  between  two  points  on  a chain. 


6.1 


PNTCHD  S 

POLYGN  S 

PRINT  S 

RESID  S 

ROSCAL  S 

STLINK  S 

UNPACK  S 

VERTEX  S 

VISIBL  S 

WHEX  S 

XCOORD  S 

XPROFL  S 

YCOORD  S 

YPROFL  E 


Compute  the  max  and  min  distances  from  a chain  to  a point. 
Determine  a polygonal  approximation  of  a chain  to  a 
specified  precision. 

Convert  chain  to  binary  plane  representation  and  print 
binary  plane. 

Compute  chain  residue. 

Rotate  and  scale  a chain. 

Store  a link  in  a packed  format  buffer. 

Unpack  a chain. 

Convert  a chain  into  lists  of  vertex  coordinates. 

Append  visible-chain-follows  signal  code. 

Compute  chain  width,  height,  U5  and  135  degree  extents. 
Append  an  x coordinate  to  a chain. 

Compute  chain  profile  along  x axis. 

Append  a y coordinate  to  a chain. 

Compute  chain  profile  along  y axis. 


The  following  labelled  common  blocks  are  defined  in  CHAP: 

CHAPMC  XYCCMP  PCHAIN 
UCHAIN  ICHAIN  WCHAIN 
PLANES  IOBUFF  STATUS 


Caution:  Note  that  users  must  not  define  these  global  names  and  common  block 
names  in  their  programs. 
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7.  DESCRIPTIONS  OF  CHAP  SUBPROGRAMS 


The  following  is  an  alphabetically  ordered  description  of  all  CHAP 
subprograms.  The  FORTRAN  data  types  of  parameters  passed  to  CHAP  subroutines 
and  functions  are  abbreviated  in  column  Type  a3  follows: 

I Integer 

R Real 

L Logical 

H Hollerith  character  string 


The  parameters  passed  to  CHAP  subroutines  are  also  classified  in 
column  headed  I/O  as  follows: 


I Input — value  is  assigned  to  this  parameter  in 

the  user's  program.  It  is  not  changed 
by  CHAP  programs . 

0 Output — value  is  assigned  to  this  parameter  in 

a CHAP  program. 

T Temporary — value  is  undetermined  before  and  after 

calls  to  CHAP  programs.  The  parameter 
is  used  for  temporary  storage  by  CHAP. 


In  the  column  headed  Dimension . a indicates  the  parameter  is  not 
an  array  but  is  a one-word  variable. 
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NAME : ANGLE 


TITLE:  Compute  angle  between  the  positive  x axis  and  line  defined  by 

two  chain  vertices. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  ANGLE ( IDENT ,NP1 ,NP2 ,A ) 


PARAMETERS: 

Name  I/O 

Type 

Dimension 

Description 

IDENT  I 

I 

- 

Chain  identification  number 

NPI  I 

I 

- 

First  chain  vertex 

NP2  I 

I 

- 

Second  chain  vertex 

A 0 

R 

Computed  angle  (in  degrees) 

COMMON  BLOCKS:  XYCOMP,  UCKAIN 
SUBROUTINES  CALLED:  UNPACK 


DESCRIPTION:  The  angle  (in  degrees)  between  the  positive  x-axis  and  the 
straight  line  defined  by  the  two  chain  vertices  NPi  and  NP2 
is  computed. 

WORKING  BUFFER  USED:  1 


NAME : APLANE 


TITLE:  Convert  and  append  a chain  into  a binary  plane. 

TYPE:  Entry  to  subroutine  BPLANE. 

CALL  FORM:  CALL  APLANE ( IDENT ) 

PARAMETERS: 

Name  I/O  Type  Dimension  Description 

IDENT  II  - Chain  identification  number 

COMMON  BLOCKS:  CHAPMC,  XYCOMP,  UCHAIN,  PLANES 

SUBROUTINES  CALLED:  Implementation-dependent.  See  appropriate  appendix. 

UNPACK 

DESCRIPTION:  This  subrout  ine  converts  a chain  into  a binary  plane  and 
stores  it  into  array  PLANE  in  the  commom  block  /PLANES/. 
The  previous  contents  of  PLANE  are  not  affected. 

WORKING  BUFFER  USED:  1 


6/77 


NAME : ARRAY 


TITLE:  Convert  chain  into  array  containing  the  links  one  per  word. 


TYPE:  Subroutine. 


CALL  FORM:  CALL  ARRAY (IDENT, LIST, N,L) 


PARAMETERS: 
Name  I/O 

TVPe 

Dimension 

Description 

IDENT 

I 

I 

- 

Chain  identification  number 

LIST 

I 

I 

N 

An  array  for  chain  links 

N 

I 

I 

- 

Dimension  of  LIST 

L 

0 

I 

- 

Number  of  chain  links  stored  in  LIST 

COMMON 

BLOCKS: 

UCHAIN 

SUBROUTINES  CALLED:  UNPACK 

DESCRIPTION:  This  subroutine  stores  the  links  of  a chain  into  array  LIST. 

One  link,  valued  0 to  7>  is  stored  into  one  word  of  the  array. 
All  signal  codes  are  skipped.  If  there  are  more  than  N links 
in  the  chain,  the  remaining  links  are  skipped  and  a warning 
message  is  printed.  L is  set  to  number  of  links  stored  in 
LIST. 

WORKING  BUFFER  USED:  3 
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NAME : AUTO 


TITLE:  Compute  the  autocorrelation  of  a chain. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  AUTO (IDENT, ACORR, PERIOD, LIST ,H) 


PARAMETERS: 

Descripticn 

Chain  identification  number 
Autocorrelation  value 
Period  of  maximum  autocorrelation 
Array  to  store  chain  directional  links 
Dimension  of  LIST 

COMMON  BLOCKS:  none 
SUBROUTINES  CALLED:  ARRAY 


Name 

I/O 

Type 

Dimension 

I 

IDENT 

I 

I 

ACORR 

0 

R 

PERIOD 

0 

I 

LIST 

T 

I 

N 

N 

I 

I 

DESCRIPTION:  This  subroutine  computes  the  maximum  autocorrelation  value 
of  a chain.  The  autocorrelation  value  ACORR  is  in  the 
range  < 0.0,1.0>.  PERIOD  is  the  number  of  links  of  the 
period  of  the  autocorrelation  function  when  ACORR  is  maximum 

WORKING  BUFFER  USED:  3 

LIMITATIONS:  Note  that  the  array  LIST  must  be  large  enough  to  store  all 
directional  links  of  the  chain  in  one  link-per-vord  format 
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NAME : BPLANE 


TITLE:  Clear  binary  plane  and  convert  a chain  into  it. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  BPLANE ( IDENT) 

PARAMETERS: 

Name  I/O  Type  Dimension  Description 

IDENT  II  - Chain  identification  number 

COMMON  BLOCKS:  CHAPMC,  XYCCMP , UCHAIN,  PLANES 

SUBROUTINES  CALLED:  Implementation-dependent.  See  appropriate  appendix. 

UNPACK 

DESCRIPTION:  This  subroutine  clears  the  contents  of  array  PLANE  in  the 
common  block  /PLANES/  and  then  converts  chain  IDENT  into 
a binary  plane  representation  and  stores  it  in  PLANE. 

WORKING  BUFFER  USED:  1 
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NAME : 


B PR I NT 


TITLE:  Print  binary  plane. 


TYPE:  Subroutine  entry  point  in  PRINT. 

CALL  FORM:  CALL  SPRINT (IW) 


PARAMETERS : 

Lane  I/O  Type  Dimension  Description 


IW  I I 


u 


Specified  window  to  be  printed 


COMMON  BLOCKS:  CHAPMC , IOBUFF,  PLANES 


SUBROUTINES  CALLED:  Implementation-dependent. 


See  appropriate  appendix. 


BPLANE 


DESCRIPTION:  A window  specified  by  IW  of  the  CHAP's  binary  plane  PLANE 

is  common  block  /PLANES/  is  printed. 

The  contents  of  IW  are 

IW(l)  - initial  x point 
IW(2)  - final  x point 
IW(3)  - initial  y point 
IW(4)  - final  y point 
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NAME : BUILD 


TITLE:  Start  to  build  a chain  by  initializing  buffer  and  inserting 

chain's  identification  number. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  BUILD ( IBUFF , IDENT ) 


PARAMETERS: 

Name  I/O 

Tyne  Dimension 

Description 

IBUFF  I 

I 

CHAP  buffer  to  store 

the  chain 

IDENT  I 

I 

Chain  identification 

number 

NAME : CENTRD 


, 

TITLE:  Compute  the  x and  y coordinates  of  the  centroid  of  a chain. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  CENTRD(IDENT,CX,CY) 


PARAMETERS 

Name 

I/O 

Type  Dimension 

Description 

IDENT 

I 

I 

Chain  identification  number 

CX 

0 

R 

Centroid  x coordinate 

CY 

0 

R 

Centroid  y coordinate 

COMMON  BLOCKS:  none 


SUBROUTINES  CALLED:  ECAREA,  MONO. 

LIMITATIONS:  Chain  IDEUT  must  be  a closed  chain. 

DESCRIPTION:  The  x and  y coordinates  (CX,CY)  of  the  centroid  of  a chain  are 
computed  as  the  ratios  of  the  chain's  moments  about  the  x and 
y axes  to  the  enclosed  area. 

WORKING  BUFFER  USED:  1 
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NAME : CHAPBD 


TITLE: 

TYPE: 

CALL  FORM: 
DESCRIPTION 


Allocate  and  initialize  the  CHAP  labelled  common  blocks. 


Block  data  program. 


Implementation  dependent.  See  appropriate  appendix. 

This  program  allocates  and  initializes  the  CHAP  labelled  common 
blocks.  These  blocks  contain  machine-dependent  constants,  chain 
buffers,  an  input/output  buffer,  binary  plane  buffer,  catalogs 
of  chains  in  memory  and  parameters  referring  to  the  sizes  of 
the  above  buffers.  Changing  these  parameters  allows  a user  to 
alter  the  sizes  of  buffers  in  the  entire  CHAP  system. 


The  following  are  common  blocks  allocated  and  initialized  by 
this  program: 

a)  Machine  dependent  constants  - block  CHAPMC 


/CHAPMC/  :iBW  tILW  .BLANK , EX  .DIGITS  ( 8 ) 


N3W 

NLW 

BLANK 

EX 


number  of  bits  per  machine  word 

number  of  octal  links  per  machi  le  word 

left  Justified  code  of  the  space  character 

left  Justified  code  of  the  character  used 
to  print  chains  from  binary  plane 


digits! 8) 

b)  Components  of  chain  directional  links  - block  XYCOMP 


left  Justified  codes  of  the  eight  octal 
digits 


/XYCOMP/  AX(8) ,AY(8) 

AX(3)  x components  of  chain  links  0 to  7 

AY(8)  y components  of  chain  links  0 to  7 
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c)  Catalog  of  chains  in  CHAP  - block  ICKAIN 

/ICHAIN/  UCH ( 9 ) * MAXNCH,MCH,CHDIR(3*MAXNCH) 

UCH(9)  not  used,  replaced  by  /STATUS/ 

MAXNCH  maximum  allowable  number  of  chains 

in  CHAP  memory 

NCH  current  number  of  chains  in  CHAP 

memory 

CHDIR( 3»MAXNCH)  catalog  of  packed  chains  in  memory: 

CHDIR(3i-2)  identification  number  of 
the  i-th  chain  in  /PCHAIN/ 

CHDIR(3i-l)  index  of  first  location  of 
the  i-th  chain  in  /PCHAIN/ 

CHDIR(3i)  index  of  last  location  of  the 
i-th  chain  in  /PCHAIN/ 

d ) Packed  chains  and  working  buffers  - blocks  PCHAIN  and  UCHAIN 
/PCHAIN/  BUFFER (PSIZE) 

BUFFER ( PSIZE ) array  for  storage  of  packed  chains 
/UCHAIN/  UNBUFF (US I ZE) 

UNBUFF(USIZE)  array  for  storage  of  working  buffers 
Both  arrays  are  to  be  dimensioned  by  the  user  according 
to  his  requirements. 

e ) Starting  indices  of  the  three  working  buffers  - clock  /CHAIN 
/WCHAIN/  WBASE(U) 

WBASE(U)  starting  indices  of  the  three  partitioned 

working  buffers  in  /UCHAIN/  UNBUFF(USIZE) . 
WBASE(k)  contains  the  dimension  of  UNBUFF 
plus  one. 
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REFERENCES: 


f ) Binary  plane  - block  PLANES 
/PLANES/  NX , NY , PLANE ( NX /NBW#NY ) 

NX  number  of  bits  (points)  in  the  x direction 

NY  number  of  bits  (points)  in  the  y direction 

PLANE  binary  plane 

g)  Input /output  buffer  - block  IOBUFF 
/IOBUFF/  IO(IOSIZE) 

IO(IOSIZE)  input /output  buffer  - usually  dimensioned 

to  132 

h)  Status  of  unpacked  chains  - block  STATUS 

/STATUS/  ID(3) ,SX(3) ,SY(3) ,SSTART(3) ,SST0P(3) ,LLINK(3) , 
PART( 3 ) ,SEOCH( 3 ) 


The  description  of  the  following  arrays  applies  to  the  chain  currently 
contained  in  each  working  buffer. 


ID  ( 3 ) 

SX(3) 

SY(3) 

SSTART(3) 

SSTOP( 3 ) 

LLINK(3) 

PART ( 3 ) 

SE0CH(3) 


the  chain  identification  number 
the  initial  x coordinate 
the  initial  y coordinate 

the  first  position  used  in  the  working  buffer 

the  last  position  used  in  the  working  buffer 

the  total  number  of  links  in  the  chain 

the  number  of  calls  to  the  routine  UNPACK  in 
the  current  unpacking  sequence  (in  the  case  of 
a chain  which  is  larger  than  the  working  buffer) 

a logical  flag  which  indicates  if  the  end-of- 
chain  indicator  was  found 


See  appropriate  appendix  for  9 listing  of  this  program. 
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NAME : CHLINE 


TITLE:  Generate  a sequence  of  links  approximating  a straight  line. 


TYPE:  Subroutine. 


CALL  FORM:  CALL  CHLINE(DELX,DELY,I ) 


PARAMETERS: 

Name  I/O 

DELX  I 

DELY  I 

I I 


Type  Dimension  Description 
I Final  x coordinate 

I Final  y coordinate 

I Specified  starting  location  in  the  array  UNBUFF 

to  store  the  link  sequence. 


COMMON  BLOCKS:  /UCHAIN/UNBUFF(USIZE) 


SUBROUTINES  CALLED:  none 


DESCRIPTION:  This  subroutine  generates  a sequence  of  chain  links,  vhich 
approximate  a straight  line  from  grid  point (0,0)  to  the 
point  (DELX, DELY  ).  The  sequence  of  links  is  stored  in  array  UNBUFF 
of  common  block  /UCHAIN/ ; locations  I through  I + MAX( |DEXL| , |DELY| )-l. 
Note  that  end-of-chain  code  is  not  appended. 

REFERENCES:  Bresenham,  J.  E.,  "Algorithm  for  computer  control  of  a digital 
plotter",  IBM  Systems  J. , (L),  1965,  pp.  25-30. 
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NAME : CHLIST 


TITLE:  List  all  chains  in  CHAP  catalog. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  CHLIST 

PARAMETERS:  none 

COMMON  BLOCKS:  CHAPMC ,PCHAIN,ICHAIII,IOBUFF 
SUBROUTINES  CALLED:  GTLINK 

DESCRIPTION:  This  subroutine  lists  on  FORTFAJI  logical  I/O  unit  6 
contents  of  the  catalog  of  chains  and  the  contents 
chains  in  memory  at  the  time  of  the  listing. 


the 

of  all 
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NAME : CHPAX 


TITLE:  Convert  a CHAP  chain  to  a PAX  chain. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  CHPAX ( ID ENT , LIST ,N ) 


PARAMETERS 

Name 

I/O 

Type 

Dimension 

DescriDtion 

IDENT 

I 

I 

- 

Chain  identification  number 

LIST 

0 

I 

N+l 

PAX  chain  array 

N 

0 

I 

— 

Number  of  chain  direction  codes  in  LIST 

COMMON  BLOCKS:  UCHAIN 

SUBROUTINES  CALLED:  UNPACK 

DESCRIPTION:  This  subroutine  converts  a CHAP  chain  into  a PAX  formatted 

list  of  directional  codes.  Such  a list  can  be  stored  into  a 
PAX  plane  by  the  PAX  routine  CHAIN.  The  PAX  direction  codes 
are  as  follows: 

6 7 8 

5 « 1 

h 3 2 

The  codes  are  stored  one  per  word  into  array  LIST  and  terminated 

by  a zero  code.  N+l  is  the  number  of  codes  in  LIST  including 
the  terminating  zero  code. 

'WORKING  3UFFER  USED:  1 


6/77 


NAME : CHPLOT 


TITLE:  Plot  a chain  on  a digital  plotter. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  CHPLOT (IDENT, STRING, SIZE) 
PARAMETERS: 


Name 

m 

Tyne 

Dimension 

Description 

IDENT 

I 

I 

Chain  file 

STRING 

I 

H 

LO  char 

Output  identification  string 

SIZE 

I 

R 

Grid  size  in  inches 

COMMON  BLOCKS:  XYCCMP.UCHAIN 

SUBROUTINES  CALLED:  Implementation  dependent.  See  appropriate  appendix. 

UNPACK 

DESCRIPTION:  This  subroutine  plots  a chain  number  IDENT  on  a digital  plotter. 

The  plotter  output  is  identified  by  a LO  character  Hollerith  string 
presumably  containing  the  user's  name  and  account  number. 

WORKING  BUFFER  USED:  1 
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NAME : CHPLOT 


TITLE:  Plot  a chain  on  a digital  plotter. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  CHPLOT  (IDENT,  GSIZE) 

PARAMETERS: 

Name  I/O  Type  Dimension 

IDENT  II 
GSIZE  I R 

COMMON  BLOCKS:  XYCOMP,  UCHAIN 

SUBROUTINES  CALLED:  UNPACK 

WORKING  BUFFER  USED:  1 

DESCRIPTION:  This  subroutine  plots  a chain  on  a digital  plotter.  The 

chain  identification  number  is  printed  at  the  start  of  the 
chain.  The  calcomp  routines  PLOTS,  PLOT,  SYMBOL,  and 
NUMBER  are  used. 


Description 

Chain  identification  number 
Grid  size  in  inches 
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NAME : COARSE 


TITLE:  Requantize  chain  on  a coarser  grid. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  COARSE  (J,XCOOR,YCOOR,N, LINKS 

PARAMETERS: 


Name 

I/O 

Type 

Dimension 

Description 

J 

i 

I 

- 

Number  of  links  in  chain 

XCOOR 

i 

T 

J+l 

X coordinates  of  chain  vertices 

YCOOR 

i 

I 

J+l 

Y coordinates  of  chain  vertices 

N 

i 

I 

- 

Number  of  times  to  expand  chain  grid 

LINKS 

0 

I 

J 

Links  of  compressed  chain 

M 

0 

I 

- 

Number  of  links  in  LINKS 

o 

o 

T* 

O 

22 

BLOCKS: 

none 

SUBROUTINES  CALLED:  none 

DESCRIPTION:  This  subroutine  compresses  a chain  by  superimposing  it  on 

a chain  grid  expanded  N times  and  then  requantizing  it.  The 
directional  links  of  the  new  chain  are  stored  one  link  per 
word  in  array  LINKS . 


NAME : CROSS 


TITLE:  Compute  the  crosscorrelation  of  two  chains. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  CROSS ( IDENT1 , IDENT2 , CORK, J, LIST ,N) 


PARAMETERS: 

Marne  I/O 

Type 

Dimension 

Description 

IDENT1 

I 

I 

% 

Chain  1 identification  number 

IDENT2 

I 

I 

- 

Chain  2 identification  number 

CORR 

0 

R 

- 

Correlation  between  chain  1 and  2 

J 

I 

I 

- 

Chain  shift  value 

LIST 

T 

T 

N 

Array  to  store  links  of  chains  1 and  2 

N 

I 

I 

- 

Dimension  of  LIST 

COMMON  BLOC 

KS: 

none 

SUBROUTINES 

CALLED: 

ARRAY 

DESCRIPTION:  This  subroutine  determines  the  degree  of  shape  similarity 

(crosscorrelation)  between  two  chains  by  comparing  the  directional 
links  of  chain  1 and  chain  2 in  pairs . The  increment 
value  J determines  the  degree  of  shift  to  the  right  of  chain  1 'with 
respect  to  chain  2.  The  result  is  the  average  pairwise  difference 
between  chains  1 and  2 for  a given  value  of  J.  The  value  of 
CORR  is  in  the  range  -1.0, 1.0  . 

LIMITATIONS:  Note  that  the  array  LIST  must  be  large  enough  to  store  all  directional 
links  of  both  chains  in  cne-link-per-word  format. 

WORKING  BUFFER  USED:  3 
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NAME : DELETE 


TITLE: 

Delete  a chain  from  CHAP 

TYPE: 

Subroutine . 

CALL  FORM: 

CALL  DELETE (IDENT) 

PARAMETERS: 

Name 

m 

Tyne  Dimension 

Description 

IDENT 

I 

I 

Chain  identification  number 

COMMON 

BLOCKS: 

ICHAIN,  PCHAIN 

SUBROUTINES  CALLED:  none 

DESCRIPTION:  This  subroutine  deletes  the  chain  number  IDENT  from  the 
catalog  of  chains  in  CHAP. 
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NAME : ECAREA 


TITLE:  Compute  the  enclosed  area  of  a closed  chain. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  EC  AREA  ( IDEIIT  ,S ) 

Description 

Chain  identification  number 
Enclosed  area 


PARAMETERS: 

Name  I/O 


I DENT 
S 


I 

0 


type  Dimension 
I 
R 


COMMON  BLOCKS:  XYCOMP,  UCHAIN 
SUBROUTINES  CALLED:  UNPACK 

LIMITATIONS:  Chain  IDENT  must  be  closed. 

DESCRIPTION:  This  subroutine  computes  the  enclosed  area  of  a chain. 


It  is  as 

sumed 

that 

the  chai 

n is  closed.  The  area 

is 

positive 

if 

it 

is 

enclosed 

in  the 

clockwise  sense 

and 

negative 

if 

it 

is 

enclosed 

in  the 

ccntraclockwise 

sense 

WORKING  BUFFER  USED:  1 
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E:  EOi 


DESCRIPTION:  The  end-of-chain  control  code  1 OUOO'  is  appended  to  chain 

IDE1TT  which  is  currently  being  built  by  CHAP,  by  placing  the 
number  10i o in  the  next  position  of  the  working  buffer. 
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NAME  : 


GTLI'IK 


TITLE:  Unpack  a link  from  a packed  format  chain. 

TYPE:  Integer  Function. 


CALL  FORM:  Variable=GTLINK( BUFFER, I) 


PARAMETERS: 

Name  I/O  Type 

BUFFER  I I 

I II 


Dimension 

M 


Description 

One  dimensional  array 

Sequential  position  of  the  digit  in 
array  BUFFER 


COMMON  BLOCKS:  CHAPMC 


SUBROUTINES  CALLED:  Implementation  dependent.  See  appropriate  appendix. 

DESCRIPTION:  This  integer  function  returns  the  T-th  digit  from  the  one 

dimensional  array  containing  links  in  packed  format , with  NLW 
octal  digits  per  word. 

LIMITATIONS:  In  order  to  stay  within  the  array  3UFFER: 

1 < I < NLU#M 
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NAME : INITXY 


TITLE:  Find  the  x and  y coordinates  of  the  chain  initium. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  INITXY (IDENT.X.Y) 


PARAMETERS : 

Name  I/O 

Tyne  Dimension 

Description 

IDENT 

I 

I 

Chain  identification  number 

X 

0 

I 

The  initial  x coordinate 

Y 

0 

I 

The  initial  y coordinate 

COMMON 

BLOCKS: 

STATUS 

SUBROUTINES  CALLED:  UNPACK 

DESCRIPTION:  This  subroutine  searches  a chain  for  the  signal  codes  that 
give  the  x and  y coordinates  of  the  initium.  Their  values 
are  returned  as  X and  Y.  If  none  are  found,  the  values  of 
X and  Y are  set  to  zero.  If  there  is  more  than  one  x or  y 
coordinate  in  a chain,  the  value  of  the  first  is  returned 
and  all  others  are  ignored. 

WORKING  BUFFER  USED:  1 


6/77 


NAME : INPUT 


TITLE:  Input  a chain  data  file. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  INPUT (ICODE) 

PARAMETERS: 

Name  I/O  Tyne  Dimension  Description 
ICODE  II  - Catalog  flag 

COMMON  BLOCKS:  CHAPMC ,PCHAIN,ICHAIH,IOBUFF 

SUBROUTINES  CALLED:  STLINK 

LIMITATIONS:  Input  chains  must  be  in  the  form  specified  below. 

DESCRIPTION:  This  subroutine  reads  one  or  more  chains  from  a data  file  which  has 
been  assigned  the  FORTRAN  logical  I/O  unit  number  5-  If  the 
parameter  ICODE  is  zero  the  catalog  of  chains  in  CHAP 
is  initialized,  else  the  chains  read  are  added  to  those 
chains  already  catalogued  by  CHAP. 


The  chain  data  files  must  be  entered  in  the  following  format : 
Columns  1-5  Chain  identification  number 

6 Ignored 

7-72  Chain  data 

73-80  Sequence  number 


In  the  chain  data  field  all  characters  other  than  octal  digits 
0 to  7 are  ignored.  This  allows  users  to  freely  insert  blanks 
or  other  characters  in  order  to  improve  the  readability  of 
chains.  The  data  file  is  terminated  by  a zero  chain  identifi- 
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cation  number. 


NAME : INTERS 


TITLE:  Find  the  intersections  of  two  chains. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  INTERS(IDENT1,IDENT2,ITYPE,INLCHA,INLCHB, NUMINT, X,Y,LEXIST,D) 

PARAMETERS : 

Name  I/O  Type 

IDENT1  I I 

I DENT 2 I I 

ITYPE  0 I 

INLCHA  0 I 

INLCHB  0 I 

NUMINT  0 I 

X T I 

Y T I 

LEXI3T  0 T 

D II 

COMMON  BLOCKS:  none 

SUBROUTINES  CALLED:  VERTEX 

LIMITATIONS:  J must  be  at  least  the  number  of  intersections.  The  dimension  D of 
arrays  X and  Y must  be  large  enough  to  store  the  coordinates  of  all 
the  links  of  both  chains,  which  is  the  number  of  links  in  both  chains 
plus  two. 

DESCRIPTION:  This  subroutine  finds  all  intersections  between  chains  1 and  2.  An 

intersection  may  be  of  two  types  - nodal  ITYPE* 1 and  non-nodal  ITYPE?2. 
The  number  of  intersections  found  is  NUMINT.  X,Y  and  LEXIST  are 
arrays  used  by  this  subroutine  to  store  the  x and  y coordinates  of 
vertices  of  the  two  chains  and  pointers  to  intersections.  INLCHA(J) 
and  INLCHB(J)  contain  the  number  of  directional  links  of  chain  1 and  2 
involved  in  the  J-th  intersection  found. 

WORKING  BUFFER  USED:  1 


Dimension  Description 

* Chain  1 identification  number 

Chain  2 identification  number 
J Type  of  intersection 

J Link  of  chain  1 in  intersection 

J Link  of  chain  2 in  intersection 

Number  of  intersections 

D Array  for  x coordinates  of  chains  1 and  2 

D Array  for  y coordinates  of  chains  1 and  2 

J Array  for  pointers  to  intersections 

Dimensions  of  X,  Y 
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NAME:  INVIS 


TITLE:  Append  invisible-chain-follows  signal  code. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  INVIS (IDENT) 

PARAMETERS: 

Name  I/O  Type  Dimension  Description 

IDENT  II  - Chain  identification  number 

COMMON  BLOCKS:  STATUS,  UCHAIN 

SUBROUTINES  CALLED:  none 

DESCRIPTION:  The  invisible-chain-follows  control  code  'OLOl'  is  appended 
to  chain  IDENT  which  is  currently  being  built  by  CHAP,  by 
placing  the  number  llio  in  the  next  position  of  the  working 
buffer . 
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NAME : LENGTH 


TITLE:  Compute  the  length  of  a chain. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  LENGTH (IDENT,CHL,LCF) 


PARAMETERS: 

Name  I/O  Tyne  Dimension 

Description 

IDENT 

II 

Chain  identification  number 

CHL 

OR 

Computed  length  of  a chain 

LCF 

II 

LENGTH  computation  flag 

COWON 

BLOCKS:  UCHAIN 

SUBROUTINES  CALLED:  UNPACK 

DESCRIPTION:  This  subroutine  computes  the  length  of  a chain.  All  even 

valued  links  are  of  length  1 unit;  all  odd  valued  links  are 
of  length  /2  units. If  the  length  of  both  visible  and  invisible 
links  is  desired  LCF  should  be  set  to  0.  For  the  length  of 
only  visible  links  LCF  should  equal  1,  and  for  only  invisible 
links  LCF  should  equal  2. 

WORKING  BUFFER  USED:  1 
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NAME : LINK 


TITLE:  Append  a link  to  a chain. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  LINK (IDENT, LINK, NA) 


PARAMETERS: 

Name  I/O 

Tyne  Dimension 

Description 

IDENT 

I 

T — 

i. 

Chain  identification  number 

LINK 

I 

I 

Link  to  be  appended 

NA 

I 

I 

Number  of  times  the  link  is  to  be  appended 

COMMON  BLOCKS:  STATUS,  UCHAIN 


SUBROUTINES  CALLED:  none 

DESCRIPTION:  This  subroutine  appends  the  directional  link  LIFK  to  chain 
IDENT  currently  being  built  by  CHAP.  The  link  LINK  is 
appended  NA  times,  where  NA  > 0. 
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NAME : LINKSQ 


TITLE:  Append  a sequence  of  links  to  a chain. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  LINKSQ( IDENT .LIST ,NL,NA ) 
PARAMETERS: 


Name 

I/O 

Type 

Dimension 

Descrintion 

IDENT 

I 

I 

- 

Chain  identification  number 

LIST 

I 

I 

NL 

An  array  of  links  to  be  appended 

NL 

I 

I 

- 

Number  of  links  in  LIST 

NA 

I 

I 

- 

Number  of  times  the  sequence  of 

COMMON 

BLOCKS: 

STATUS 

.UCHAIN 

in  LIST  is  to  be  appended 

SUBROUTINES  CALLED:  none 

DESCRIPTION:  This  subroutine  appends  the  sequence  of  ML  directional 

links  stored  in  array  LIST,  one  link  per  word,  to  chain  IDENT 
currently  being  built  by  CHAP.  The  sequence  of  links  is 
appended  HA  times,  where  NA  > 0. 
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NAME : LM0M1 


TITLE:  Compute  first  moment  of  a chain  about  a line. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  LMOMl( IDENT ,X1  ,Y1  ,X2 ,Y2 ,FMNT) 


PARAMETERS: 

Name  I/O 

Type 

Dimension 

Description 

I 

I 

- 

Chain  identification  number 

XI 

I 

R 

- 

X coordinate  of  first  point 

of 

line 

Y1 

I 

R 

- 

Y coordinate  of  first  point 

of 

line 

X2 

I 

R 

- 

X coordinate  of  second  point 

of 

line 

Y2 

I 

R 

- 

Y coordinate  of  second  point 

of 

line 

FNCJT 

0 

P. 

- 

Computed  first  moment 

COMMON 

BLOCKS: 

none 

SUBROUTINES  CALLED:  LNAP.EA , CENTRD 
LIMITATIONS:  Chain  IDENT  must  be  a closed  chain. 

DESCRIPTION:  This  subroutine  computes  the  first  moment  of  a chain  about 
a line  defined  by  points  ( XI, Y1  ) and  ( X2,Y2  ). 

WORKING  BUFFER  USED:  1 
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NAME : LM0M2 


TITLE:  Compute  the  second  moment  of  a chain  about  a line. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  LM0M2( IDENT, XI, Y1,X2,Y2,SMNT) 


PARAMETERS 

Name 

IDENT 

I/O 

Type 

i 

Dimension 

Description 

Chain  identification  number 

XI 

I 

R 

- 

X coordinate  of 

first  point 

Y1 

I 

R 

- 

Y coordinate  of 

first  point 

X2 

I 

R 

- 

X coordinate  of 

second  point 

Y2 

I 

R 

- 

Y coordinate  of 

second  point 

SMUT 

0 

R 

- 

Computed  second 

moment 

COMMON  BLOCKS:  XYCOMP.UCHAIN 
SUBROUTINES  CALLED:  INITXY, UNPACK ,MOM2A 

DESCRIPTION:  This  subrouti  ne  computes  second  moment  of  a chain  about 
a line  defined  by  points  (Xl,Yl)  and(X2,Y2). 

LIMITATIONS:  Chain  IDENT  must  be  a closed  chain. 

WORKING  3UFFER  USED:  1 
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NAME : LNAREA 


TITLE:  Compute  the  area  between  a chain  and  a line. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  LNAREA(IDENT,X1,Y1,X2,Y2,AREA) 


PARAMETERS: 
Name  I /O 

IDENT  I 

Tyne  Dimension 

I 

Descrintion 

Chain  identification  number 

XI 

I 

R 

X coordinate  of 

first  point 

Y1 

I 

R 

Y coordinate  of 

first  point 

X2 

I 

R 

X coordinate  of 

second  point 

Y2 

I 

R 

Y coordinate  cf 

second  point 

AREA 

0 

R 

Computed  area 

COMMON 

BLOCKS: 

XYCOMP ,UCHAIN 

SUBROUTINES  CALLED:  ECAREA.INITXY, UNPACK, MOMIA.CENTRD 
LIMITATIONS:  Chain  IDEtJT  must  be  a closed  chain. 

DESCRIPTION:  This  subroutine  computes  the  area  between  a chain  and  a 
line  defined  by  points  (X1,Y1)  and  (X2,Y2)  . 

WORKING  BUFFER  USED:  1 
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NAME : LNCHD 


TITLE:  Compute  the  max  and  min  distance  from  a chain  to  a line. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  LHCHD(IDENT,X1,Y1,X2,Y2,CLDMaX,CLDMIH,JMAX,JIC[N) 


PARAMETERS 

Name 

I/O 

Type 

Dimension 

Description 

IDENT 

I 

I 

- 

Chain  identification  number 

XI 

I 

R 

- 

X coordinate  of  first  line  point 

Yl 

I 

R 

- 

Y coordinate  of  first  line  point 

X2 

I 

R 

- 

X coordinate  of  second  line  point 

Y2 

I 

R 

- 

Y coordinate  of  second  line  point 

CLDMAX 

0 

R 

- 

Maximum  distance  between  chain  and 

CLDMIN 

0 

R 

- 

Minimum  distance  between  chain  and 

JMAX 

0 

T 

1 

- 

Number  of  link  at  maximum  distance 

JMIN 

0 

I 

— 

Number  of  link  at  minimum  distance 

COMMON  BLOCKS: 

none 

SUBROUTINES  CALLED:  INITXY,  ARRAY 

DESCRIPTION:  This  subroutine  finds  the  minimum  and  maximum  distances  between 
a chain  and  a line  defined  by  points  (XI, Yl)  and(X2,Y2). 

J?1AX  and  JMIU  are  numbers  of  directional  links  at  the  maximum 
and  minimum  distances,  respectively. 

WORKING  BUFFER  USED:  1,3. 
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NAME : MATCH 


TITLE:  Match  two  chain  profiles. 

TYPE:  Subroutine. 


CALL  FORM: 

CALL 

MATCH ( XX , YY , PR0F1 

,N1,PR0F2,N2) 

PARAMETERS 

Name 

I/O  ; 

Type 

Dimension 

Description 

XX 

0 

I 

- 

Relative  position  in  the  x 

direction 

YY 

0 

I 

- 

Relative  position  in  the  y 

direction 

PR0F1 

I 

I 

2,N1 

First  profile  sequence 

PR0F2 

I 

I 

2 ,N2 

Second  profile  sequence 

N1 

I 

I 

- 

Dimension  of  PR0F1 

N2 

I 

I 

- 

Dimension  of  PR0F2 

' 

i 

COMMON  BLOCKS : none 


SUBROUTINES  CALLED:  none 

DESCRIPTION:  This  subroutine  computes  the  best  match  of  two  profile 

sequences  PRGFL1  and  PR0FL2.  XX  and  YY  are  the  number  of 
grid  positions  to  shift  the  first  chain  in  the  X and  Y 
direction  for  the  optimum  match. 

REFERENCES:  Freeman,  H. , "On  the  tempi ate -layout  problem",  N.Y.U. 

Technical  Report  CRL-39,  February  1975. 
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N AME : MAXMIN 

TITLE:  Compute  the  max  and  min  extents  of  a chain  in  x and  y directions. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  MA»HN(IDENT,XY) 


PARAMETERS: 

Name  I/O 

Tyne 

Dimension 

Description 

IDENT  I 

I 

- 

Chain 

XY  0 

I 

1* 

Array  for  the  min/max  values 

COMMOM  BLOCKS:  XYCOMP.UCHAIN 
SUBROUTINES  CALLED:  UNPACK 


DESCRIPTION:  This  subroutine  computes  the  maximum  and  minimum  coordinates 
of  chain  number  I DENT  in  the  x and  y directions.  The 
values  are  stored  into  the  array  XY(L)  in  the  following 
order: 

XY(l)  - maximum  x coordinate  of  chain 
XY(2)  - minimum  x coordinate  of  chain 
XY(3)  - maximum  y coordinate  of  chain 
XY(L)  - minimum  y coordinate  of  chain 


WORKING  BUFFER  USED:  1 


NAME : M0M1 


TITLE:  Compute  the  first  moment  of  a chain  about  0,^5, 90,135  axes. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  MOM! ( IDENT , DEGREE , FMNT ) 


PARAMETERS: 

Description 

Chain  identification  number 
Axis  in  degrees 
Computed  first  moment 

COMMON  BLOCKS:  XXCOMP,UCHAIN 


Name  I , C 
IDENT  I 
DEGREE  I 
FMNT  0 


Type 

I 

I 

R 


Dimension 


SUBROUTINES  CALLED:  INITXY, UNPACK 


DESCRIPTION:  This  subroutine  computes  the  first  moment  of  a chain 
about  0,  U5,  90,  and  135  degree  axes.  The  axis  is 
selected  by  setting  DEGREE  to  0,  L5,  90,  or  135- 

LIMITATIONS:  The  chain  IDENT  must  be  a closed  chain. 

WORKING  BUFFER  USED:  1 


, 


1 
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NAME : M0M2 


TITLE:  Compute  the  second  moment  of  a chain  about  0,1+5,90,135  axes. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  MQM2 ( IDENT , DEGREE ,SMNT ) 


PARAMETERS 

Name 

I/O 

Tyne  Dimension 

Description 

IDENT 

I 

I 

Chain  identification  number 

DEGREE 

T 

T — 

Selected  axis  in  degrees 

SMUT 

0 

R 

Computed  second  moment 

COMMON  BLOCKS: 

XYCOMP,  UCHAIN 

SUBROUTINES  CALLED:  INITXY,  UNPACK 

DESCRIPTION:  The  second  moment  of  a chain  is  found  as  the  product  of 

the  area  of  the  chain  and  the  square  of  the  distance  from 
the  axis.  Thus  the  x and  y components  of  each  element  with 
respect  to  the  given  axis  is  found.  The  moment  of  an  area 
encircled  in  a clockwise  sense  above  the  x axis  or  in  a 
counterclockwise  sense  below  the  axis  is  positive. 

LIMITATIONS:  Chain  IDENT  must  be  a closed  chain. 

WORKING  BUFFER  USED:  1 
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NAME : If  PACK 


TITLE:  Pack  a chain. 


TYPE:  Subroutine. 


CALL  FORM:  CALL  NPACK ( IDENT , NB , START , STOP ) 


PARAMETERS 

: 

Name 

I/O 

Type 

Dimension  Description 

IDENT 

i 

I 

Chain  identification  number 

NB 

i 

I 

Number  of  working  buffer 

START 

i 

I 

Index  of  the  first  location  of  the 
unpacked  chain 

STOP 

i 

T 

Index  of  the  last  location  of  the 
unpacked  chain 

COMMON  BLOCKS:  CHAPMC,ICHAIN,PCHAIN,UCHAIN 
SUBROUTINES  CALLED:  STLINK.GTLINK 

DESCRIPTION:  This  subroutine  creates  a new  packed  chain  from  the 

contents  of  working  buffer  NB  which  contains  a chain  in  the 
unpacked  format.  START  is  the  index  of  the  first  used 
location  in  the  working  buffer  NB.  STOP  is  the  last 
location  used  in  the  working  buffer  NB.  IDENT  is  the 
identification  number  of  the  new  chain  . 
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NAME : OCTNUM 


TITLE:  Form  an  octal  number  from  a sequence  of  links. 

TYPE:  Integer  function. 

CALL  FORM:  Variable=OCTIIUM( BUFFER, I, N) 


PARAMETERS: 

Name  I/O 

£ge. 

Dimension 

DescriDtion 

BUFFER  I 

i 

M 

One-dimensional  array 

I I 

i 

- 

Number  of  the  first 

digit 

H I 

i 

- 

Number  of  digits  to 

form  the  number 

COMMON  BLOCKS: 

CHAPMC 

SUBROUTINES  CALLED:  Implementation  dependent.  See  appropriate  appendix. 

DESCRIPTION:  This  function  returns  the  integer  value  (decimal)  of  a number 
formed  from  N octal  digits  stored  in  one  dimensional  array 
BUFFER  as  described  for  programs  GTLINK  and  STLINK.  The  I-th 
octal  digit  is  in  the  most  significant  position  and  the 
(l+N-l)-th  octal  digit  is  in  the  least  significant  position  of  the 
number . 

LIMITATIONS:  In  order  to  stay  vithin  the  array  BUFFER: 

lil  - NLW#M 
1-N  - NLW»M-I+1 
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NAME : OUTPUT 


TITLE:  Output  chains  to  I/O  unit  7 in  format  read  by  INPUT. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  OUTPUT ( IDENT ) 


PARAMETERS: 

Name  I/O 

IDENT  I 


Type  Dimension  Description 
I - Chain  identification  number 


COMMON  BLOCKS:  CHAPMC ,PCHAIN,ICHAlN,IOBUFF 


SUBROUTINES  CALLED:  GTLINK 


DESCRIPTION: 


This  subroutine  writes  a chain  on  a data  file,  which  has 
been  assigned  FORTRAN  logicaj.  I/O  unit  number  7.  The 
chain  data  file  will  be  output  in  the  following  format: 


Columns 

1-5 

Chain  identification  number 

6 

blank 

7-72 

Chain  data 

73-80 

Sequence  number 

The  sequence 

numbers  begin  at  lOOio  and  are 

incremented  by  lOOio- 


6/77 


NAME : PD I ST 


TITLE:  Compute  the  distance  between  two  points  on  a chain. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  PDIST(IDENT,P1,?2,DIST) 

PARAMETERS: 

Name  I/O  Type 

IDENT  I I 

PI  II 

P2  II 

DIST  0 R 

COMMON  BLOCKS:  XYCOMP.UCHAIN 
SUBROUTINES  CALLED:  UNPACK 

DESCRIPTION:  This  subroutine  computes  the  straight  line  distance  between 
two  points  on  a chain  defined  by  the  vertices  PI  and  P2. 

WORKING  BUFFER  USED:  1 
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NAME : PNTCHD 


TITLE:  Compute  the  max  and  min  distances  from  a chain  to  a point. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  PNTCHD ( ID ENT , XP , YP , DMAX , DMIN , LMAX , LMIN ) 

PARAMETERS: 

Name  I/O  Type  Dimension 

IDENT  II 

XP  I R 

YP  I R 

DMAX  OR 

DMIN  OR 

LMAX  0 1 

LMIN  01 

COMMON  BLOCKS:  XYCOMP.UCHAIN 
SUBROUTINES  CALLED:  INITXY, UNPACK 

DESCRIPTION:  This  subroutine  computes  the  maximum  and  minimum  distances 
between  a chain  and  a point  defined  by  (XF,YP).  LMAX  and 
LMIN  are  sequential  numbers  of  chain  directional  links  at 
maximum  distance  DMAX  and  at  minimum  distance  DMIN, 
respectively. 

WORKING  BUFFER  USED:  1 


Description 

Chain  identification  number 
X coordinate  of  a point 
Y coordinate  of  a point 
Computed  maximum  distance 
Computed  minimum  distance 
Number  of  link  at  maximum  distance 
Number  of  link  at  minimum  distance 
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NAME : POLY  ON 


TITLE:  Determine  a polygonal  approximation  of  a chain  to  a 

specified  precision. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  POLYGN  (IDEir, XCOORD, YCOORD, ICL, IOP, TOL, IC) 


PARAMETERS 

Name 

I /0  Type 

Dimension 

Description 

IDENT 

I I 

- 

Chain  identification  number 

XCOORD 

0 I 

J 

X coordinates  of  chain  vertices 

YCOORD 

0 I 

J 

Y coordinates  of  chain  vertices 

ICL 

0 I 

J 

Pointers  to  vertices  of  polygon 

IOP 

I I 

J 

Stack  used  in  approximation 

TOL 

I R 

- 

Tolerance  of  the  polygonal  approximation 

IC 

0 I 

- 

Number  of  vertices  of  polygon 

COMMON  BLOCKS : XYCOMP 

,UCHAIN 

SUBROUTINES  CALLED:  INITXY, UNPACK 

DESCRIPTION:  This  subroutine  computes  the  polygonal  approximation  of  a 

a chain.  XCOORD  and  YCOORD  store  the  coordinates  of  the  chain 
vertices.  ICL  stores  pointers  to  the  vertices  of  the  polygon,  i.e., 
XCOORD(ICL(J) ) and  YCCGED(ICL(J) ) are  the  coordinates  of  the  J-th 
vertex  of  the  polygon,  IC  is  the  number  of  vertices  of  the  polygon. 
TOL  is  set  in  a calling  program  to  the  required  tolerance  for  the 
polygonal  approximation. 

LIMITATIONS:  J is  the  number  of  vertex  points  in  the  chain. 

WORKING  BUFFER  USED:  1 
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REFERENCES: 


Ramer,  U. , "An  Iterative  Procedure  for  the  Polygonal 
Approximation  of  Plane  Curves",  Computer  Graphics  and 
Image  Processing,  1,  1972,  pp.  2kk-256 . 
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NAME : PRINT 


TITLE:  Convert  a chain  to  binary  plane  representation  and  print 

the  binary  plane. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  PRINT ( IDENT,  IW) 


PARAMETERS: 

Name  I /0 

Type 

Dimension 

Description 

IDENT  I 

I 

- 

Chain  identification  number 

IW  I 

T 

1 

h 

Window  to  be  printed 

COMMON  BLOCKS:  CHAPMC ,IOBUFF, PLANES 

SUBROUTINES  CALLED:  Implementation-dependent . See  appropriate  appendix. 

DESCRIPTION:  This  subroutine  converts  the  chain  IDENT  into  a Dinary  plane 
and  then  outputs  a window  of  the  plane  as  specified  by  IW 
to  FORTRAN  logical  unit  number  6. 
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NAME : RESID 


TITLE:  Compute  chain  residue. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  RESID ( IDENT ,LRES1 ,NLRES1 ,LRES2 ,NLRES2 ) 


PARAMETERS: 
Name  I/O 

IDENT  I 

LRES1  0 

NLRES1  0 

LRES2  0 

NLRES2  0 


Type  Dimension 
I 
I 
I 
I 
I 


Description 

Chain  identification  number 
Number  of  first  type  links  in  residue 
Number  of  links  LI  in  residue 
Number  of  second  type  links  in  residue 
Number  of  links  L2  in  residue 


COMMON  BLOCKS:  XYCOMP.UCHAIN 


SUBROUTINES  CALLED:  UNPACK 


DESCRIPTION:  This  subroutine  computes  the  links  in  the  residue  of  a 
chain.  There  are  at  most  two  types  of  links  in  the 
residue,  which  differ  by  unity,  i.e.  LRE31=LRES2+l(mcd  8). 
Each  type  of  residue  link  occurs  0 or  greater  number  of 
times  as  specified  by  NLRES1  for  type  LRES1  and  NLRES2 
for  type  LRES2  of  residue  links. 


WORKING  BUFFER  USED:  1 
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NAME : ROSCAL 


TITLE:  Rotate  and  scale  a chain. 

TYPE:  Subroutine 


CALL  FORM:  CALL  ROSCAL  ( IDENT1 , IDENT2 , ANGLE , XSCALE,  YSCALE) 


PARAMETERS 

• 

Name 

I/O 

Type  Dimension 

Description 

IDENT1 

T 

I — 

Identification  number  of  the 
be  rotated  and  scaled 

chain 

to 

IDENT2 

I 

I 

Identification  number  of  the 
chain 

result 

ing 

ANGLE 

I 

R 

Rotation  angle  in  degrees 

XSCALE 

I 

R 

Scale  factor  in  x direction 

YSCALE 

I 

R 

Scale  factor  in  y direction 

COMMON  BLOCKS: 

UCHAIN,  V CHAIN 

SUBROUTINES  CALLED:  UNPACK,  NPACK,  INITXY,  3UILD,  LINK 


WORKING  BUFFER  USED:  1,3 

DESCRIPTION:  This  subroutine  rotates  and  scales  a chain.  The  resulting 
chain  is  stored  with  the  specified  identification  number. 

If  the  angle  is  a multiple  of  90°  and  the  scale  factors  are 
both  1.0,  then  a very  efficient  algorithm  is  used. 
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NAME : STL INK 


TITLE:  Store  a link  in  a packed  format  buffer. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  3TLINK( LINK, BUFFER, I) 

PARAMETERS: 

Name  I/O  Type  Dimension  Description 
LINK  II  - Octal  digit 

BUFFER  II  M One- dimensional  array 

I II  - Desired  sequential  position  of  LINK  in 

BUFFER 

COMMON  BLOCKS:  CHAPMC 

SUBROUTINES  CALLED:  Implementation  dependent.  See  appropriate  appendix. 

DESCRIPTION:  This  subroutine  inserts  the  octal  digit  LINK  into  the 
one  dimensional  array  BUFFER  as  the  I-th  sequential  octal 
digit  in  packed  format  vith  NLW  octal  digits  per  word. 

LIMITATIONS:  In  order  to  stay  within  the  array  3UFFER: 

1-1  - NLW#M 
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NAME : UNPACK 


P 


TITLE:  Unpack  a chain. 

TYPE:  Subroutine. 


CALL  FORM:  CALL  UNPACK (IDENT.NB, FIRST, ISTART,ISTOP,EOCH) 


PARAMETERS 

Name 

I/O 

Type 

Dimension 

DescriDtion 

IDENT 

x 

x 

- 

Chain  identification  number 

NB 

I 

I 

- 

Number  of  working  buffer 

FIRST 

I 

L 

- 

First  part  of  chain  to  be  unpacked 

I START 

0 

I 

- 

Index  of  the  first  location  of  the 
unpacked  chain 

ISTOP 

0 

I 

- 

Index  of  the  last  location  of  the  1 
chain 

EOCH 

0 

L 

_ 

End  of  chain  flag 

COMMON  BLOCKS:  CKAPMC,PCHAIN,UCHAIN,ICHAIN,WCHAIN 


SUBROUTINES  CALLED:  QTLINK ,OCTNUM 

DESCRIPTION:  This  subroutine  unpacks  a portion  of  a packed  chain 
into  working  buffer  NB.  I3TART  is  the  index 
of  the  first  used  location  in  the  working  buffer  N3,  ISTOP 
is  the  index  of  the  last  used  location  in  the  working  buffer 
NB.  The  logical  flag  EOCH  is. TRUE. if  the  end  of  chain 
indicator  was  reached. 

The  subroutine  ?s  called  multiple  times  to  ur.nack  a long  chain. 

FIRST  is  to  be  set  to  .TRUE,  in  the  call  if  the  first  part  of  a chain 
is  to  be  unpacked.  It  is  to  be  set  to  .FALSE,  in  the  call  to  unpack 
additional  portions  of  the  chain  until  EOCH  is  .TRUE. 
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NAME : VERTEX 


TITLE:  Convert  a chain  into  lists  of  vertex  coordinates. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  VERTEX ( IDENT , XCOORD , YCOORD ,N ,L ) 


PARAMETERS 

Name 

I/O 

Type 

Dimension 

Description 

IDENT 

I 

I 

- 

Chain  identification  number 

XCOORD 

0 

I 

N 

An  array  of  x coordinates 

YCOORD 

0 

I 

N 

An  array  of  y coordinates 

N 

I 

I 

- 

Dimension  of  XCOORD  and  YCOORD 

L 

0 

I 

- 

Number  of  vertex  coordinates  in 
XCOORD  and  YCOORD 

COMMON  BLOCKS:  XYCOMP ,UCHAIN 

SUBROUTINES  CALLED:  INITXY, UNPACK 

DESCRIPTION:  This  subroutine  converts  a chain  into  two  arrays.  Array 
XCOORD  contains  the  x coordinates  of  the  vertex  points  of 
the  chain  and  array  YCOORD  contains  the  y coordinates. 

Note  that  XCOORD(l)  and  YCOCRD(l)  contain  the  initial 
coordinates  of  the  chain.  If  there  are  more  than  IT— 1 
directional  links  in  the  chain  the  remaining  links  are 
skipped  and  a warning  message  is  printed.  L is  set  to  the 
number  of  coordinate  pairs  in  the  XCOORD  and  YCOORD  arrays. 

FORKING  BUFFER  USED:  1 
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NAME : VT3IBL 


TITLE:  Append  visible-chain-follows  signal  code. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  VISIBL(IDENT) 

PARAMETERS : 

Name  I/O  Type  Dimension  Description 

ID ENT  II  - Chain  identification  number 

COMMON  BLOCKS:  STATUS .UCHAIN 

SUBROUTINES  CALLED:  none 

DESCRIPTION:  The  visible-chain-follows  signal  code  '0L02'  is  appended 
to  chain  IDENT  which  is  currently  being  built  by  CHAP, 
by  placing  the  number  12io  in  the  next  location  of  the 
working  buffer. 
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NAME : WHEX 


TITLE:  Compute  chain  width,  height,  1*5  and  135  degree  extents. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  WHEX(IDENT,ITYPE,W) 

PARAMETERS: 

Name  I/O  Type  Dimension  Descrintion 
IDENT  II  - Chain  identification  number 

ITYPE  II  - Type  of  computation  desired 

W OR  - Computed  value 

COMMON  BLOCKS:  UCHAIN 
SUBROUTINES  CALLED:  UNPACK 

DESCRIPTION:  This  subroutine  computes  the  width,  height,  1+5  degree 
extent  or  135  degree  extent  of  a chain  by  specifying 
ITYPE  to  be 

1 - width 

2 - height 

3-U5  degree  extent 
1*  - 135  degree  extent 

WORKING  BUFFER  USED:  1 


6/77 


NAME : XCOORD 


TITLE:  Append  an  x coordinate  to  a chain. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  XCOORD( IDENT, X) 


PARAMETERS 

Name 

m 

Type 

Dimension  Description 

IDENT 

i 

I 

Chain  identification  number 

X 

i 

I 

X coordinate 

COMMON  BLOCKS:  STATUS, UCHAIN 
SUBROUTINES  CALLED:  none 

DESCRIPTION:  The  x coordinate  value  and  signal  code  is  appended  to 

a chain  with  identification  number  IDENT  which  is  built 
by  CHAP. 

LIMITATIONS:  0-X  - 32768! 0 
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NAME : XPROFL 


TITLE:  Compute  a chain  profile  along  x axis. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  XPROFL ( IDENT , FLAG, PROFL, II, LIST ,NN) 


PARAMETERS: 


Name 

I/O 

Dimension 

Description 

IDENT 

T 

1 

T 

1 

- 

Chain  identification  number 

FLAG 

L 

- 

Flag  indicating  upper  profile  (.TRUE.) 
or  lower  profile  (.FALSE.)  desired 

PROFL 

0 

I 

2,N 

Array  of  profile  pairs 

N 

I 

I 

- 

Dimension  of  array  PROFL 

LIST 

T 

j 

NN 

Array  to  store  chain 

NN 

I 

I 

_ 

Dimension  of  array  LIST 

COMMON  BLOCKS : XYCOMP 

SUBROUTINES  CALLEO:  ARRAY, MAXMIN 

DESCRIPTION:  This  subroutine  forms  profile  pairs  of  a chain-coded  template. 

The  profile  pairs  are  computed  relative  to  the  enclosing 
rectangle  made  of  x and  y maximum  and  minimum  positions  of  the 
chain.  When  FLAG  is  true  the  profile  is  computed  along  the 
upper  edge,  when  false  the  profile  is  computed  along  the  lower 
edge  of  the  template. 

REFERENCES:  Freeman,  H. , "On  the  template-layout  problem",  N.Y.U. 

Technical  Report  CRL-39,  February  1975 

WORKING  BUFFER  USED:  1,3 
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NAME : YCOORD 


TITLE:  Append  a y coordinate  to  a chain. 

TYPE:  Subroutine. 

CALL  FORM:  CALL  YCOORD ( IDENT , Y ) 


PARAMETERS 

Name 

I/O 

Type  Dimension 

Description 

IDENT 

I 

T 

JL 

Chain  identification  number 

Y 

I 

I 

Y coordinate 

COMMON  BLOCKS:  STATUS, UCHAIN 
SUBROUTINES  CALLED:  none 

DESCRIPTION:  The  y coordinate  and  signal  code  is  appended  to  a chain 
with  identification  number  IDENT  which  is  currently  being 
built  by  CHAP. 

LIMITATIONS:  0-Y  - 3276810 
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8.  AN  EXAMPLE  OF  A CHAP  USER  PROGRAM 


A simple  CHAP  user  program  is  listed  on  the  following  page.  The 
program  reads  an  input  chain  data  file  which  consists  of  2 records.  The 
first  contains  a '1'  in  column  5 and,  starting  in  column  7,  the  chain 
"121230007651*55il0U00" . The  second  record  consists  of  a zero  in  column  5. 

The  expected  output  is  shown  in  the  program  listing.  The  chain  is  listed, 
then  the  image  of  the  chain  is  printed  from  a 1*1x6  point  window  of  the  CHAP 
binary  array  as  specified  by  window  size  IW.  Then  the  enclosed  area,  length, 
first  moment  about  the  y axis  and  second  moment  about  the  x axis  of  the 
chain  are  computed  and  printed.  The  output  of  the  example  program  is 
given  on  page  8.3. 
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V. 

AN  EXAMPLE  OF  CHAP  PROGRAMMING 


I MPL I C I T INTEGER ( A-Z ) 

REAL  L f S f FMNT  * SMNT 
DIMENSION  I LK 4 ) 

DATA  IW  / 0>40>0>5  / > LCF  / 0 / 

READ  IN  ONE  SAMPLE  CHAIN y WITH  I DENT . NO.  = 1 
CALL  INPUT (0) 

ID  = 1 

LIST  THE  CHAIN 
CALL  CHL 1ST 

• 

PRINT  THE  CHAIN 
CALL  PRINT ( I D r IW) 

FIND  THE  ENCLOSED  AREA 
CALL  ECAR'EA  ( ID » S ) 

WRITE <6 , 100)  IDrS 

100  FORMAT ( 31H-THE  ENCLOSED  AREA  OF  CHAIN  N0.»I3»4H  IS  ,F7.2/) 

FIND  THE  LENGTH  OF  THE  CHAIN 
CALL  LENGTH ( ID , L r LCF ) 

WRI TE<  6 v 1 10 ) ID  i>  L 

110  FORMAT ( 24H  THE  LENGTH  OF  CHAIN  NO.,  13,  4H  IS  r F7.2/) 

FIND  THE  FIRST  MOMENT  ABOUT  THE  Y AXIS 
CALL  MOM 1 ( ID  * 90 » FMNT ) 

WRITE < 6 ? 120  > FMNT 

120  FORMAT < 31 H THE  FIRST  MOMENT  ABOUT  Y AXIS: > F7.2/> 

FIND  THE  SECOND  MOMENT  ABOUT  THE  X AXIS 
CALL  M0M2  < ID  * 0 * SMNT ) 

WRITE (6, 130)  SMNT 

130  FORMAT < 32H  THE  SECOND  MOMENT  ABOUT  X AXIS : r F7.2/) 


>C 

C 

>C 

C 

>C 


EXPECTED  OUTPUT 
CHAIN  LISTING- 
ENCLOSED  AREA- 
LENGTH- 

FIRST  MOMENT  < Y 
SECOND  MOMENT- 


12123  00076  54504  04  00 

11.50 

17.90 

AXIS)  -32.17 
121.08 


STOP 

FNn 
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ft EXECUTION  PEG I NS 


CONTENTS  OF  THE  CHAIN  BUFFER- 


CHAIN  ID  NUMBER  = 1 

12123  00076  54554  0400 


END  OF  BUFFER 


CHAIN  FILE  NO.  1 

XXXX 
X X 
X X 
X XX 

XX 


THE  ENCLOSED  AREA  OF  CHAIN  NO.  1 IS  11. 
THE  LENGTH  OF  CHAIN  NO.  1 IS  17.90 
THE  FIRST  MOMENT  ABOUT  Y AXIS:  -32.17 
THE  SECOND  MOMENT  ABOUT  X AXIS.'  121.06 

■ftEXECUTION  TERMINATED 

ft 
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10.  APPENDICES 


A.  Multiple-Computer  Version 

B.  I3M  360  Version 


Appendix  A.  MULTIPLE-COMPUTER  VERSION 


Most  of  the  CHAP  subroutines  are  totally  independent  of  the  computer 
on  vhich  CHAP  is  to  be  run.  All  that  is  required  is  a standard  ANSI 
FORTRAN  IV  compiler.  The  only  exceptions  are  the  chain  packing  and  unpack- 
ing routines,  vhich  depend  on  the  word  length  of  the  computer,  as  well  as 
some  minor  variations  in  the  particular  FORTRAN  versions.  It  is  planned 
to  maintain  a single  master  version  of  CHAP.  This  version  will  consist 
of  the  'onion  of  all  the  statements  required  for  all  the  different  computer 
systems  on  vhich  CHAP  is  to  be  implemented.  It  will  not  be  an  executable 
version.  Special  control  cards,  identified  by  the  letter  "M"  in  column  1 
vill  be  used  to  indicate  the  computer  systems  to  which  the  statements  that 
follow  apply.  For  each  computer  system,  a unique  column,  beginning  with 
column  10,  is  assigned  in  the  M-cards.  An  entry  (any  alphanumeric 
character  other  than  'blank')  in  one  or  more  columns  indicates  that  all  the 
cards  that  follow  (until  the  next  M-card  is  encountered)  apply  to  all  the 
computer  systems  that  have  an  entry  in  their  assigned  column.  (The  entries 
may  be  coded  to  facilitate  readability). 

Statements  that  apply  to  all  computer  system  versions  are  preceded  by 
an  M-card  that  either  contains  an  entry  in  every  assigned  computer  column 
beginning  with  column  10,  or  contains  an  A in  column  9-  In  the  latter  case, 
any  entries  in  columns  10  and  above  are  ignored. 

To  generate  a CHAP  version  for  a particular  computer  system,  the 
master  version  is  first  pre-processed  by  a small  program  to  (l)  delete  all 
statements  for  which  the  preceding  M-card  either  does  not  have  an  A in 
column  9 or  does  not  have  an  entry  in  the  column  corresponding  to  the 


The  result  will  be 


desired  computer  version,  and  (2)  delete  all  M-cards. 
an  executable  FORTRAIJ  IV  program  that  will  run  on  the  desired  computer 


system. 


The  master  version  will  be  centrally  maintained,  revised,  and 
extended  as  found  necessary  or  desirable. 


Appendix  B.  IBM  360/370  VERSION 


This  version  of  CHAP  was  written  for  an  1314  360  computer  and 
compiled  by  a FORTRAN  (G)  compiler.  It  was  tested  on  the  IBM-360/67 
at  RPI,  using  the  Michigan  Terminal  System. 

Three  IBM  assembly  language  integer  functions  are  used.  They  are 
LAND,  IOR,  and  SHIFT,  and  are  described  on  the  following  pages.  They  are 
used  in  the  process  of  packing  and  unpacking  chains.  The  memory 
organization  of  CHAP  packed  chains  is  shown  in  the  fig-are  below.  Ten  (10) 
links  are  packed  in  one  32-bit  word  in  this  version.  Bits  0 and  1 of  each 
word  are  not  used. 


BITS /WORD  ( NBW ) 

LINKS /WORD  (NLW) 

32 

10 

0-1.2.  31 

0-1  2 

31  0-1  2 

i - l • 1 T 

i LjLi j . L 

■m 

1 T 1 T 1 

1.  L l L__  i 

: 1 1 10 

HBI 

• 20  21 i 

i i i 

W0RD1  W0RD2  WORD 3 


The  Block  Data  program  CHAPBD  must  be  link-edited  with  the  user's 
main  program  to  allocate  the  required  memory  space  correctly.  The  work- 
ing buffers  (unpacked)  are  allocated  600  words  in  total,  and  the  packed 
suffer  is  allocate!  500  words  in  this  version.  The  I3M  version  of  CHAPBD 
is  listed  in  this  appendix. 
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NAME : YPROFL 


TITLE:  Compute  chain  profile  along  y axis. 

TYPE:  Entry  to  subroutine  XPROFL. 


CALL  FORM: 

: CALL  YPROFL (IDENT, FLAG 

, PROFL, N, LIST, NN) 

PARAMETERS: 

Name  I/O 

Type 

Dimension 

Description 

ID  ENT 

I 

I 

- 

Chain  identification  number 

FLAG 

I 

L 

- 

Flag  indicating  right  (.TRUE.)  or  left 
profile  (.FALSE.)  is  desired 

PROFL 

0 

I 

2,N 

Array  for  profile  pairs 

N 

I 

I 

- 

Dimension  of  array  PROFL 

LIST 

TO 

I 

NN 

Array  to  store  chain 

NN 

I 

I 

- 

Dimension  of  array  LIST 

COMMON  BLOCKS:  XYCOMP 


SUBROUTINES  CALLED:  ARRAY , MU 

DESCRIPTION:  This  subroutine  forms  profile  pairs  of  a chain-coded  template. 

The  profile  pairs  are  computed  relative  to  the  enclosing 
rectangle  made  of  x and  y maximum  and  minimum  positions  of 
the  chain.  When  FLAG  is  true  the  profile  is  computed  along 
the  right  edge,  when  false  the  profile  is  computed  along  the 
left  edge  of  the  template. 

REFERENCES:  Freeman,  H. , "On  the  template-layout  problem",  N.Y.U. 

Technical  Report  CRL-39,  February  1975 

WORKING  BUFFER  USED:  1,3 
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ADDITIONAL  LIST  OF  GLOBAL  NAMES 


The  following  is  a list  of  the  global  names  which  are  defined  in 


the  IBM-360/370  version  of  CHAP.  They  are  ail  integer  functions 

and  are 

described 

on  pages  B.3-B.5- 

Name 

Type 

Description 

LAND 

IF 

The  logical  and  of  two  variables. 

IOR 

IF 

The  logical  or  of  two  variables. 

SHIFT 

IF 

Left  circular  shift  of  an  integer  variable. 

The  following  list  is  of  the  CHAP  routines  which  call  the 
implementation  dependent  routines  and  the  routines  they  call. 

CHAP  routine  Routines  called 


A PLANE 

SHIFT, 

IOR 

3PLANE 

SHIFT, 

IOR 

3PRINT 

CKPLOT 

SHIFT, 

I AND 

SHIFT, 

LAND 

GTLINK 

OCTNUM 

SHIFT , 

IAND 

PRINT 

SHIFT, 

LAND 

STLINK 

SHIFT, 

IAND 

3.2 


NAME : IAIID 


TITLE:  Form  bit-by-bit  logical  AND  of  two  variables. 

TYPE:  Integer  function. 


CALL  FORM:  Variable  = IANT>(WORDl,  WORD2) 


PARAMETERS: 
Name  I/O 

WORD1  I 

WORD2  I 


Tyne  Dimension  Description 

I - An  integer  variable 

I - An  integer  variable 


COMMON  BLOCKS : none 

SUBROUTINES  CALLED:  none 

DESCRIPTION:  IBM  assembly  language  routine  for  the  I3M  360 

implementation  of  CHAP.  Assembled  by  ASSEMBLER  (F). 

The  logical  AND  of  the  two  parameters  is  returned. 
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NAME : IOR 


TITLE:  Fora  bit-by-bit  logical  OR  of  two  variables. 

TYPE:  Integer  function. 

CALL  FORM:  Variable  = I0R(W0RD1,  W0RD2 ) 

PARAME'I  ERS : 


Name 

I/O 

Type 

Dimension  Description 

W0RD1 

I 

I 

An  integer  variable 

W0RD2 

I 

I 

An  integer  variable 

COMMON 

3L0CKS : 

none 

SUBROUTINES  CALLED:  none 

DESCRIPTION:  IBM  assembly  language  routine  for  the  IBM  360 

implementation  of  CHAP.  Assembled  by  ASSEMBLER  ( 

The  logical  OR  of  the  two  parameters  is  returned. 


B.U 


6/77 


NAME : 


SHIFT 


TITLE: 


TYPE: 


CALL  FORM: 


Left  circular  shift  of  an  integer  variable. 

Integer  function. 

Variable  = SHIFT  (WORD,N) 


PARAMETERS : 

Name 

I/O 

Type  Dimension 

WORD 

I 

I 

N 

I 

I 

Description 

An  integer  variable 

Number  of  bits  to  shift  0<N<30 


COMMON  3L0CKS : none 

SUBROUTINES  CALLED:  none 

DESCRIPTION:  Left  circular  shift  of  the  30  lower  bits  of  an  integer 

variable  is  performed.  The  two  leading  bits  of  the  result 
are  set  to  zero. 

IBM  assembly  language  routine  for  the  IBM  360 
implementation  of  CHAP.  Assembled  by  ASSEMBLER  IF). 
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3.5 


oCOCK  data  ChAPbD 


CHAP  COMMON  bLOc'K  ALLOCATION  AND  INITIALIZATION 

THIS  BLOCK  DA  I A FRGCNAM  ALLOCATES  AND  IMTI AL  I ZES 
ALL  COMMON  bLULKS  CF  CHAP 

THIS  CLOCK  UrfU  PROGRAM  MUST  BE*  COMPILEC  OR 
LIMn-IOI  rm;  1 iigtT  HLK  * I 1 H Tut:  MAIN  PfiGGhAV 


IMPLICIT  INTCGEB(A-Z> 


COMMON  /CHAPN'C/  M3W.  NLA  , 0 X . 6L  AN  K , D I C.  I TS 
COMMON  /XYCJMP/  AX,  AY 
COMMON  /PCHAIN/  L'UFLER 

COMMON  /LCHmIN/  ONEUFh 

COMMON  /(CHAIN/  uch.maxnch.nch, choir 
COMMON  /-.iCHAIN/  «BASc 
COMMON  /PLANTS/  NX , NY • PLANT 

COMMON  /I09QFF/  IQ 

C U M M On  / S r A I ub/  IU.SX.SY.sST'ART",  SSTUP ,1  L I NK , PAP  I , Sf 

CHA«AC1ER  CODES  CF  THE  CCTAL  DIGITS  C TC  7 

Cl  MENS  IJN  O I C II  b ( '!  ) 

X AND  y COMPONENTS  fF  CHAIN  I.  1 NKS 
OIMENjIJN  AX(H) »AY(B) 

CHAIN  BUFFERS 
DIMENSION  BUFFER  ( ‘3  00  ) 

DIMtNSIJN  Chi)  IP  ( 60  ) 

WORKING  U 0 F F F.  R S 
DIMENSION  UN0UFFI6C0J 
ElNAfiy  PLANE 

OIMENaiJN  PLANt(^HO)  ' “ 


INPOT/HUTPLT  bUFFt R 
OI MENS  ION  10(131) 

D I ME  No  I On  oCTim 
U l MENS  ION  V»6ASE(4) 


sh 


STATUS SwcS Yo 
LOGICAL  SECCH<3J 

DIMENSION  I 0 ( 3 ) ,SX I 3)  , 3Y ( 3 ) , 3 ST  ART (3)  . SSTCP ( J ) . LL I NX ( J ) , PARTI  3) 


MAChINL  uLFc'KUl  NT  CL  NS  JAMS 

NOW  NUMbE?  OF  BITS  P 

NEW  MJWpfcP  OF  LINKS 

DATA  NCvn.NLrt  /JO,  10/ 

TTTTl  TTTTTT  IN  THE  IBM  VERSION  ct-  c HAP  ETT3  C 
rtMt  USED  TO  STORE  THE  LINKS  OR  BINARY  PLANE 


NUMbE?  OF  BITS  PER  WCRD 
NUMBER  OF  L INKS  PER  WCRD 

"n  cF  cEaP  EIT5  c-J  1 


CHARACTER  symeols 

TmH  ex  ,ul/vne-,"d * c i rs  "? — nrrv  rT^crr7~,T",~,~2r' . ■ s ■ . ■ > 

BASF  INDEXES  ur  The  f HP  E E WORKING  EOFFEPS 
DATA  *v  EASE  / I . 2 J 1 .40  1 , oO  l / 

X AND  Y CCLKbl  NATES  OF  LINKS 

DATA  A A / 1,1  , G • 3 ■ I » I , U , 1 / , A Y / 0,  1 *1  , I ,ij»  — l * — * I »“T 
PICTURE  PLANE 
DATA  NX,  NY  / 120, 12  0/ 


DATA  OCH  .MAXNCH  .NCH.CHJ  IR  / y * 0 , 2 0 , 0 , 0*  0/ 

STATUS  VALLES 

DATA  I D , SX  , S Y , Sb  T ART  , ..  STOP  ill  I N<  , F ART  / 2 1 * 0 / 
UATA  SECCH  / 3<>.  TRUE./ 

END 


